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Abstract 

We present a calculation of the polarized p meson cloud in a nucleon using 
time-ordered perturbation theory in two different variants advocated in the lit- 
erature. We calculate the induced difference between the distributions Ad{x) 
and Au{x). We use a recent lattice calculation to motivate an ansatz for the 
polarized valence quark distribution of the p meson. Our calculations show that 
the two theoretical approaches give vastly different results. We conclude that 
Ad{x) — Au{x) can be of relevant size with important consequences for the com- 
bined fits of polarized distribution functions. 

It is experimentally established beyond doubt that d{x) and u{x) in the nucleon 
differ substantially [|I], 0, ^. The difference can be described as due to the meson cloud 
in the nucleon. The relevant process is shown in fig. 0. The nucleon splits up into a 
meson and a baryon and the meson is hit by the photon in a deep inelastic scattering 
experiment. 

The dominant role in the unpolarized case is played by the pion. It is clear that a 
proton prefers being accompanied by a vr"*" instead of an vr" since in the first case the 
remaining baryon can be a neutron, but in the second case the remaining baryon with 
the lowest lying mass is the A(1232) resonance. Therefore one can find an enhancement 
of d quarks over u quarks in the proton. A large number of papers has been published 
on this subject, e.g. g |, |, 0, |, ||, |ig. 

Proceeding a step further in this direction one can ask the question whether there is 
also an asymmetry in the polarized sea, more precisely between the distributions Ad{x) 
and Au{x). Up to now it is not possible to extract such an asymmetry from experiments 
but the situation could change in the future. It is clear that the meson cloud model 
predicts such an asymmetry. The mechanism is the same as in the unpolarized case, 
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Figure 1: The MNB process. A nucleon splits up into a meson and a baryon and 
the meson is hit by a photon in a deep inelastic scattering experiment. 



now with the lightest polarizable meson — the p-Meson with a mass of 770 MeV — in 
the role of the pion. There have already been attempts to derive this asymmetry from 



quark distribution functions calculated in the chiral quark-soliton model |]TT], |I2|. 

In this publication we present a calculation of the difference between the light 
polarized antiquark distributions. We take into account the processes where a nucleon 
fluctuates into p and A^ [pNN) and into p and A (pNA). Let us consider the process 
in fig. |1| in an infinite momentum frame with P, P', k and q being the four-momenta 
of the incoming nucleon, the outgoing baryon, the meson and the photon respectively. 
We are only interested in spin 1 mesons and S, S' and A shall denote the helicities of 
the nucleon, the baryon and the meson. 

The contribution 6q^ of the MNB process in fig. |l| to the distribution of antiquarks 
of a given flavour with helicity cr in a nucleon A^ is described by the well known 
convolution formula lOll 



xdq (X) 



[ dyfMNB{y)j<fM,\ 



X _„ I X 

y 



where 



X = (2) 
2P-q ^ ' 



is the usual Bjorken variable with Q = —q and 

k-q 
y=p-q 

is the momentum fraction of the nucleon carried by the meson. fMNsiv) denotes the 
probability distribution for a meson M with helicity A in the nucleon due to the MNB 
process and q^x gives the probability to find an antiquark of a given flavour with 
helicity a in a Meson with helicity A. As a shorthand notation for the convolution we 
use 6q°' = fMNB ® 'ill X following. Now let us consider a proton with hehcity 
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+ 1/2 (or briefly |) as the incoming nucleon. With the usual definitions of the helicity 
asymmetry in the nucleon 

Aq = - q^ (4) 

and in the meson 

AqM = Aglf,i - Ag[^^i = Agj^^ - Agl^__, (5) 
and with the notation 

^fuNB = fllNB - ImNB (6) 

for the polarization of the meson cloud we can write e.g. the contribution of the process 
p — s> p+ + n to the polarization of the d quarks as 

S{Ad) = Afp+pn^Adp+. (7) 

Charge symmetry tells us that 

Adp+ = Aup- = 2AdpQ = 2AupO =: Avp. (8) 

Taking into account all possible isospin channels of the pNN and pNA processes we 
obtain 

S{Ad) = (^^AfpNN + l^fpNA^ ® Avp (9) 
S{Au) = (^^AfpNN + ^A/p^A j ® Avp (10) 



and finally 



Ad-Au=[ '^AfpNN - ^AfpNA ] ® AVp. (11) 



fpNN and fpNA are the total p distributions summed over all isospin channels with a 
proton as the initial nucleon. 

For the calculation of the p distributions we use time-ordered perturbation theory 
(TOPT) in the infinite-momentum frame |jl3|, [1^ which guarantees on-shellness for the 
mesons in contrast to the usual covariant approach where the meson is not on its mass 
shell. The TOPT approach has become very popular in recent times since the fact 
that the meson is no longer off-shell removes a big problem of the early calculations. 
As a disadvantage there is no energy conserving 5-function at the MNB vertex. 

In the infinite momentum frame the three momentum of the proton is (0,0, p) 
with p tending to infinity and the momentum of the meson can then be written as 
{k ± cos ip,k^ sin ip^yp) with a transverse momentum k^. It has been shown that only 
time-ordered diagrams with all particles moving forward in the infinite momentum 
frame contribute, others are suppressed by 1 /p^. Therefore y takes only values between 
and 1 ^ iMf. One can write the meson distribution as an integral 

ftiNBiy) = J2 dkl\^rs^{y,kl)\' (12) 

s' "^0 
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where is the hehcity dependent amphtude for the MNB process considered. The 
hehcity s' of the baryon is not observed, so it is summed over. The amphtudes are 
given by @] 

The masses of the nucleon, the outgoing baryon and the p meson are denoted here 
with M, M' and m and 



^MBiy, ki) = + (14) 



''1. 

y y 



is the invariant squared mass of the meson-baryon state. The vertex function V^f/'^ 
contains the vector and spinor structure of the vertex contracted with the fields and 
depends on the Lagrangian used for the effective meson nucleon interaction. As usual 
a form factor must be introduced for the interaction which is parameterized in this 
case in exponential form 



G{y, k]_) = exp 



(15) 



as a function of the squared mass smb and with the cutoff parameter Amnb- 



The Lagrangians of the effective pNN and pNA interaction are |1T5[ 

^pNN = 9pNN'ip7''(pf,'4' + fpNN'ipcr'"'ip{df,(j)^ - d^(f)f,), (16) 

and 

J^pNA = 9pNA'ipil5Y'^''{dp(j)u - d^(f)f,) + h.c. (17) 

with the Dirac field ip for the nucleon, the Rarita-Schwinger field for the A and the 
vector field (p^ for the p meson. The vertex functions therefore are 

y's' = 9pNNu{P',s'h%{k,X)u{P,l/2) + 

+ 2fpNNu{P\ S')ia^''k^tlu{P, 1/2), (18) 
V^} = gpNAlf'iP', Sy^,Yu{p, 1/2) [tk^elik, A) - ik,e;ik, A)] , (19) 

where m is a Dirac spinor, If^ a Rarita-Schwinger spinor and e'^ the polarization vector 
of the p meson. 

There were some discussions whether this is the correct choice |T^. The reason 
for this is the appearence of some explicit factors of k in (|1^) and (0). They arise from 
the derivatives of the meson field in the interaction Lagrangians. The argument is that 
in spite of the absent energy conservation one should not use the on-shell momentum 
vector k of the meson in that case (A) but replace it by the vector k = P — P' instead 
(B). Note that only the space-like components of both vectors coincide. The definition 
of y and the polarization vectors of the meson remain unchanged. We do not want 
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Figure 2: Distributions fpN^iy) of the p mesons from the pNN process for hehcities 
A = 1 (dashed hne), (dotted hne) and —1 (dashed-dotted hne). The sohd hne 
shows the sum over all helicities and the line with circles the polarization AfpN^iy) = 
ffiNwiy) " fpNNiy)- Method (A) and a cutoff ApN]\f = 900 MeV in the formfactor are 
used. 

to make a decision in this matter. Obviously such calculations still lack a really solid 
theoretical basis and have to be interpreted as effective models to be judged by their 
phenomenological success. Both methods have been used to calculate unpolarized 
contributions of the p meson to quark distributions, see e.g. ^ for (A) and |TB[ for 
(B). In this respect both procedures do roughly equally well. 

We will present here calculations for both methods to get a feeling for the theoretical 
uncertainties involved. We find in fact large differences between the numerical results 
indicating that the theoretical treatment needs substantial improvement. Combining 
all the input we obtain for the helicity dependent and isospin summed distribution 
functions from the pNN process the results given in the appendix for both prescrip- 
tions. For the numerical calculations we used the coupling constants gpj^j^/{4:7i) = 0.84, 
fpNN = 6.lgpNN/i4:M) and g^p^Ji^Ti) = 20.45 given in [|, Figs. | and | show 
the distributions of the p mesons from the pNN and the pNA process respectively 
for different helicities obtained from method (A). In Figs. ^ and ^ the same curves are 
plotted from method (B). 

Note that method (A) favours the p mesons with helicity +1 without suppressing the 
other helicities too much. In contrast, method (B) leads to an overwhelming dominance 
of helicity +1 for the pNN channels but for pNA we have more mesons with helicity 
— 1. Note also that method (B) gives a peak at x ~ 0.55 for the unpolarized distribution 
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Figure 3: The same as in fig. H, now for tlie pNA process obtained from method (A) 
and with a cutoff Ap^A = 900 MeV. 



from the pNN process while the distributions calculated with method (A) peak below 
0.5. From the mass ratio of the p on one side and the nucleon or the A on the other 
side we would naively expect the meson to carry less than 50% of the total momentum 
in average. 

In Figs. ^ and0 we compare the polarized distributions Af{y), which are relevant 
for our problem, for both methods and different cutoffs. The cutoffs for the pNN 
and pNA vertices are not very well known. Holtmann et al. use A = 1100 MeV 
for the pNN vertex and A = 980 MeV for the the pNA vertex from fits to diffractive 
scattering data. With view to the serious problems of fixing the ttNN formfactor 
we decided to vary the cutoff parameters and keep the unpolarized p distributions of 
both methods in comparable size. The unpolarized curves can also be compared with 
those in |^ for method (A) and in for method (B). Unfortunately the dependence 
on the cutoff is for the p much larger than for the pion as already stated by Thomas 
and Melnitchouk 0. Figs. ^ and |^ show the large differences between (A) and (B). 
One notes also that both models predict a partially negative polarization of the p cloud 
from the pNA process. 

A big obstacle for considering contributions of mesons to the polarized sea of the 
nucleon is up to now our lack of knowledge about the polarized structure functions 
of vector mesons [18|. It is common to set the unpolarized valence quark distribution 
of the p equal to that of the pion. For our purpose we need the helicity asymmetry 
Avp{x) of a valence quark in the p meson which has not yet been measured. A recent 
lattice calculation ||l9l shows that the polarization of the valence quarks in the p should 
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Figure 5: The same as fig. 0, now obtained from metfiod (B) witli a cutoff parameter 
ApTVA = 850 MeV. 
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Figure 6: The polarization Afpj^]\r{y) of the p meson cloud from the pNN process 
calculated with method (A) and cutoff parameters Ap^rTv = 1100 MeV, 1000 MeV, 
900 MeV and 800 MeV (solid lines from top to bottom) and with method (B) and the 
cutoffs ApNN = 950 MeV, 850 MeV, 750 MeV and 650 MeV (dashed lines from top to 
bottom). 



be approximately 60%. Therefore we set Avp{x) = 0.6Vt^{x) where Vj^ is the valence 
quark distribution of the pion. Here we use the parameterizations of Gliick, Reya and 
Vogt for the pion [|^. This is a rather bold ansatz but it enables us to get a first 
estimate for the flavour asymmetry in the polarized sea of the nucleon. 

In figs. H we present our final results. The contributions to xAd{x) — xAu{x) 
from the pNN and pNA processes are calculated with prescriptions (A) and (B). For 
both methods we chose four different cutoff parameters A = Ap^vAf = ^pNA- The first 
point which attracts attention is the large difference between the results of (A) and 
(B). In view of the wide range of possible results it is clear that the influence of the 
polarized p mesons is not necessarily negligible. Current parameterizations for xAq{x) 



give typically a value around 0.02 or less for x = 0.1 |21|, Thus the asymmetry 
in the polarized nucleon sea could be a 50% effect. References |jTT|, |T2| predict even a 
larger asymmetry. Note that in these papers Au{x) and Ad{x) are defined with the 
opposite sign ||23|. Tab. || gives our numerical values for 



0.01 



dx {Ad - Au) (20) 



for the different models plotted in fig. p. 

Let us summarize our work. We calculated the helicity dependent distributions 



8 



0.15 




-0.05' ' ' ' ' ' 

0.2 0.4 0.6 0.8 1 

y 

Figure 7: The same as fig. § now for tlie pNA process and again cutoff parameters 
ApTVA = 1100 MeV, 1000 MeV, 900 MeV and 800 MeV (solid lines from top to bottom) 
for method (A) and ApivA = 950 MeV, 850 MeV, 750 MeV and 650 MeV (dashed lines 
from bottom to top at y = 0.4) for method (B). 

of p mesons in a proton with two models. Under convolution the different shapes of 
the p distributions do not translate directly into the quantity xAd{x) — xAu{x) but 
remarkable differences remain. Also the determination of the correct formfactors from 
unpolarized data should be improved. Bearing all that caveats in mind it is clear that 
our results show that the flavour asymmetry in the polarized sea can be relevant for 
the present DESY, SLAC and CERN spin experiments. 

We strongly thank G. Filler, W. Weise and M. Gockeler for discussions. The fruit- 
ful conversation with A. W. Thomas, W. Melnitchouk and A. Szczurek is gratefully 
acknowledged. This work was supported by BMBF. 



(A) 


(B) 


/ 


A [MeV] 


/ 


A [MeV] 


0.0064 


1100 


0.027 


950 


0.0052 


1000 


0.011 


850 


0.0030 


900 


0.0034 


750 


0.0013 


800 


0.00067 


650 



Table 1: Results for the Integral / = J^^-^ dx {Ad{x) — Au{x)) for both methods and 
different formfactors A = Ap^N = Ap^vA- 
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Figure 8: The asymmetry x{Ad{x) — Au{x)) calculated with the p distributions of 
method (A) and cutoff parameters Ap^viv = \na = HOO MeV, 1000 MeV, 900 MeV 
and 800 MeV (sohd lines from top to bottom) and with method (B) and ApjvAr = 
ApTVA = 950 MeV, 850 MeV, 750 MeV and 650 MeV (dashed lines from top to bottom) 
at g2 ^ 2 GeVl 



Appendix 



Here we give our results for the p distributions. Using the on-shell meson momentum 
k at the vertex (A) we obtain 



fpNNiy) 



fpNNiy) 



dklF\y,k] 



87rV(l-l/)^yo {M^-sl^{y,kl)f 

[g' [kl + y'M'] + +igfyM [kl + y {yHP - (1 - y)m')] + 
+ Af WM'^kl + (y^M^ - (1 - y)m^y 

3 r dklF^{y,kl) 



(21) 



167rV(l-?/)'m2 7o {M^-sl^{y,kl)f 

[kl + y^M' - (1 - y)m^] ' + Agfy'm^M [kl + y^M' 
- (1 - y)m'] + Afm' [(2 - yfkl + y'M^] }, 



(22) 
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.-1 ^ 3 dklF\y,kl) 

kl[9\i - yf - ^fgy{i - y)M + 4f [kl + y^M^] } 

for the pNN process and 

H . . ^ g' r dklF\y,kl) 

hNA[y) 247rV(l - iM^-sl^{y,kl)) 

{kl + kiM'' [3-4yil-y)] + 

+ A;iM' [Ay{l - yfm^M + 2y^M'3 + y'^{2 - yfW^] - 
- 27/2(1 - yfm^MM'^ + y^M'^ + (1 - y)6m^M2|, 

^/W r dklF\y,k] 

jpNA[y) 



+ A;i [2M'2 + (1 - yf{M^ + M'^)] + 
+ y'^M'^ - 2y\l - y)MM'^ + y''{l - yfM'^M'^y 

r dklF\y,kl) 



fpNAiy) 247rV(l-y)^M'2 7o {M^-sl^{y,kl)f 

|A;1m2 + kl ['iy^M'^M'^ - 4y(l - y)m''MM' + (1 - yfm^] + 
+ 3y*M^M'^ - 67/2(1 - y)m^MM'^ + 3(1 - yfm^M'^y 

for the pN/S. process. 

For the vector k — P — P' (B) the results are 



3 f 

fi^Ay) - 87rV(i - yY I 



dklF\y,k] 



Sn^y^l-yfJo (M^ - sl^{y,kl)f 
{g' [kl + yV] + +4gfM [y{l + y)kl + 2y'M'] + 

+ 4/2 [kl + hkly'^W + Ay^M^] }, 

rfA;iF2(y,A;i) 



fpNNiy) — 



lQnW{l-yYm^ Jo {M^-sl^{y,kl)f 



,212 



[kl + y^M^ - {l-y)rrY]' 
[/(I - y)2 - 2^/7/2(1 - y)M + f {kl{2 - yf + yV)] , 
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3 r dklF\y,kl) 

87rV(l - y)' io {M^-sl^{y,kl)f 

kl[9\i - yf - ^fgyii - y)M + Af [kl + y^M^] }, 

9^ r dklF\y,kl) 

247r2y3(i _ y)^M'^ (^2 _ ,2^^y^ ^2 ))2 

{kl + ki [(1 - y^M"^ - 4(1 - yfyMM' + 

+ (4y2 -Ay + 3)M''] + kl [2y^M'^ + y\2 - y fU'^ - 

- 2(1 - yfyM^ + 4(1 - yfy^M^M' - 2(1 - yfy^MM'^ + 
+ 2(1 - y)^yM^M'^] + (1 - y)YM^ - 2(1 - yfy^M^M''' - 

- 2(1 - y fy^M^M'^ + (1 - y fy^M'^M'^ + 

+ 2(1 - y fy^MM'^" + |/^M'<^}, 

487r2y3(l _ y)4^2^/2 _ ^2^(^^ ^2 ))2 

{fcl + fci[(|/2 + 3)M'2 + (1 - y)(l - 3y)M^ - 

- 2(1 - yW] + A;l [(((1 - y)M2 - M'')y + (1 - y)m2)' - 

- 2(((1 - y)M^ - M'^)y + (1 - y^) 

((1 - yfM^ + (2 + (1 - yf)M'^) + y^M'^ - 

- 2(1 - y)y^MM'^ + (1 - yfy'^M^M'^] + 

+ A;l [(((1 - y)M' - M'^)y + (1 - y)m'Y 
{{l-yfM^ + {2 + {l-yf)M'^)- 

- 2(((1 - y)M'^ - M'^)y + (1 - y)m^) {y^M'^ - 2(1 - y)y^MM' 
+ (1 - y)VM2M'2)] + [((1 - y)M' - M")y + (1 - y)m']' 
[y'^M'^ - 2(1 - y)y^MM'^ + (1 - ?/)^?/^M2M'2] }, 

("^ dklF\y,kl) 
247rV(l-2/)'M'2X (M2-4^(y,A;i))' 

+ A;l [4yMM' + y^M^ + (1 - y)^^^^ + 

+ 3M'2 + 2|/M'2] + A;i [(1 - yfy^M" - 4(1 - y)y^M^M' + 

+ (y^ + Qy)M'^ -2{-y^- 4y^ + 3|/)M^M'^ + 

+ lOy^MM'^] + Sy'^M"^ + 6y^MM'^ + 

+ ?,y^M'^M'^ - 6(1 - y)y^M^M'^ + 3(1 - yfy'^M^M''^ - 
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+ 6{1 -y)y^M^M"^y (32) 
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